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We review a concept of the Moscow potential (MP) of the NN interaction. On the 
basis of this concept we derive by quantum inversion optical partial potentials from 



U 

o r 

the modern partial-wave analysis (PWA) data and deuteron properties. Point-form 

C*~> ! (PF) relativistic quantum mechanics (RQM) is applied to the two-body deuteron 

CN ' photodisintegration. Calculations of the cross-section angular distributions cover 

> 

photon energies between 1.1 and 2.5 GeV. Good agreement between our theory and 



recent experimental data confirms the concept of deep attractive Moscow potential 
with forbidden S- and P-states. 

PACS numbers: 12.39.Pn, 13.40.-f, 13.75.Cs 



INTRODUCTION 



X 

Opportunities to observe manifestations of quark degrees of freedom in nuclear reactions 
at intermediate energies attract attention of scientific community for a long time. It was 
noted [l| that the most suitable subject of research here is the deuteron as the simplest 
nucleus where the secondary rescattering has little effect on the primary process. 

The deuteron photodisintegration at photon energies of ~ 2 GeV generates great interest 
among experimentalists [l[ Jy, \i\ and theoreticians 0, [f3, Q, |g[ with the main emphasis 
on the properties of the NN system which are beyond the scope of realistic mesonic NN 

n n q 

potentials [6( and can be interpreted within quark concepts [5|, [7(. First, it was shown in 
papers by Khar'kov group [6| that starting from mesonic potentials it is possible to explain 
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the g?7 — > data at energies E 7 > 1 GeV only if a revision of electromagnetic part of 
the theory is done and instead of the ordinary nucleon electromagnetic form factors the 
essentially different ones are used with poles of third order. Second, the phenomenological 



theory of Regge poles was taken as a basis in Refs. 



5| with selection of dominant poles 



according to the quark string model [h]]. Free parameters of these theories make it possible 
to describe the experimental data reasonably well. Third, also giving reasonable results 
the hard rescattering model was developed [7| within a semiempirical approach, when the 
photon is absorbed by a quark of one of the nucleons and then the hard rescattering of this 
quark by another nucleon takes place. The wave function amplitude of the final np-state 
with large relative momentum is evaluated empirically by extrapolation of the corresponding 
np-scattering experimental data. 

In this paper we use the PF RQM to treat the deuteron photodisintegration in a Poincare- 
invariant way. Modern development of the RQM and exhaustive bibliography is presented in 
the review by B.D. Keister and W. Polyzou 111]. The PF is one of the three forms proposed 
by Dirac [l2]. Other two are the front form and instant form. These forms are associated 
with the different possibilities for putting interactions in generators of the Poincare group. 
All the forms are unitary equivalent [13| but each has certain advantages. Most of the 
calculations in nuclear physics have been performed in the instant and front form. Only 
in recent years important simplifying features of the PF were realized. These features are 
connected with the fact that in the PF all the generators of the homogeneous Lorentz group 
are free of interactions. Thus only in the PF the spectator (impulse) approximation (SA) 
preserves its spectator character in any reference frame [jj], [3]. For an electromagnetic 
NN process the SA implies that the NN interaction does not affect the photon-nucleon 
interaction and therefore sum of the one-particle electromagnetic current operators may be 
taken as an electromagnetic current operator for the system of interacting nucleons. It is 
supposed that the SA may be valid when the process is quick due to the large momentum 
transfer. General covariant PF expressions for the electromagnetic current operator for 



composite systems are given in Refs. 
factors of various composite particles 



15 



18 



16]. The PF SA was applied to calculate form 



19j with reasonable results. In our calculation 



(f 

of the proton-proton bremsstrahlung 20j it was shown that the PF SA violates the continuity 
equation for the NN current operator, but the violation is relatively small for the considered 
kinematics. 
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In this paper we show that recent deuteron photodisintegration data at E 1 = 1.5—2.5 GeV 



4j confirm the Moscow NN potential model 2l| characterized by deep attractive partial 



potentials with forbidden S- and P-states. In this study the Moscow partial potentials are 



reconstructed from the iViV PWA data within the energy range < Eiai, < 3 GeV 24]. 



This reconstruction is based on our approach to the inverse-scattering problem for optical 



potentials 22]. 



The plan of the paper is as follows. In Sec. II, we review a concept of the Moscow potential 
(MP) of the NN interaction. In Sec. Ill, we present the optical Moscow-type NN potential 



derived by quantum inversion [22|] within the relativistic quasipotential approach JjJ, 123 ]. 



We show that the modern PWA data of NN scattering 



24 



are compatible with the concept 



of the MP. In Sec. IV, the formalism of PF RQM [U|, |15j, |20( is applied to the high-energy 
energy deuteron photodisintegration. Results and future prospects are discussed in Sec. V. 
In Appendix A we give necessary details of the calculation techniques. In Appendix B in the 
PF SA we derive an expression for the momentum Qn transferred to the nucleon and show 
that Qn is not the same as the momentum transfer seen by the deuteron. The expression 
is a generalization of the similar expression for the elastic electron-deuteron scattering jl^]. 



II. POTENTIALS WITH FORBIDDEN STATES IN NUCLEAR PHYSICS 

In description of systems of composite X particles consisting of some y particles it is a 
common approach to exclude explicit degrees of freedom of y particles. In simplest case of 
the XX system the microscopic Hamiltonian that includes all possible pair yy interactions 
is substituted by an effective Hamiltonian (by sum of X particle kinetic energy terms and 
of an effective XX potential). The common requirement is that the effective Hamiltonian 
would give for the XX system the same spectrum and the same corresponding relative 
motion wave functions as the initial microscopic Hamiltonian. In some cases the effective 
Hamiltonian has redundant eigenvalues and eigenstates, which must be disregarded. These 
eigenstates are called forbidden states and the effective XX potential is called then "the 
potential with forbidden states". 



For instance, in the oscillator shell model of the potential theory of a — a scattering 26] 
the antisymmetric wave function of the 8 Be nucleus ground state (eight-nucleon configura- 
tion s 4 p 4 and orbital permutation symmetry [f] x = [44]) being projected onto a — a channel 
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results in 4S-wave relative motion wave function (see our review [27|). This wave function 
accumulates all four oscillation quanta of the system and has two nodes. Momentum distri- 
butions corresponding to such wave functions were investigated in quasielastic knock-out of 
a particles from p-shell nuclei by intermediate energy photons 28J]. The OS- and 2S-states 
of a — a relative motion are forbidden as far as they correspond to the lower s 8 and s 6 p 2 
eight-nucleon configurations respectively, which are forbidden by the Pauli principle. Basing 
on these considerations, a concept of the deep attractive a — a potential with OS, 2S and 
2D forbidden bound eigenstates was elaborated j^ . According to the concept, there is 
no repulsive core in the a — a interaction and a particles can penetrate into each other. 
Forbidden bound eigenstates take lowest energy levels. Unforbidden eigenstates (including 
scattering ones) being orthogonal to the forbidden eigenstates have nodal structure at short 
range. For instance, the S-wave relative motion wave function has two nodes in the region of 
a — a overlap. This model is substantiated by the phase shift analysis based on the general- 
ized Levinson theorem (GLT) 29j. For example, the S- wave phase shift of a — a scattering 
equals 360° at zero energy, rises up to 540° at the energy slightly above the low-lying 4S- 
resonance and then runs down with increasing energy within the broad energy range up to 
Eiab ~ 200 MeV where the phase shift approaches the asymptotic region of small values and 
becomes negative due to absorption [30|. Such picture of the S-wave phase shift behavior 



was confirmed by experiments performed in a broad energy range |3l|], while .D-wave phase 
shift behavior shows one forbidden state. Phase shifts of higher waves do not show forbidden 
states (see [30|] for further details). 

In case of the NN system, the concept of the deep attractive NN potential with forbidden 
states appeared in 1975 [21] when we analyzed the pp-scattering phase shift data extended at 
the time up to = 6 GeV. It was shown that the singlet S'-wave phase shift data with an 
extended gap between low- and high-energy groups of data can be interpolated by a smooth 
curve if the empirical low-energy group is raised 180°. This interpolation demonstrates 
decrease of the S'-wave phase shift in the broad energy range from zero up to Ei a b > 5 GeV 
as a manifestation of the GLT. The high-energy part (Ei a b — 3 — 6 GeV) of the interpolation 
for the S-wave remains in the asymptotic region of small values, corresponding to the Born 
approximation. The energy dependence of the singlet D-wave phase shift is smooth and there 
is no need to raise the initial values. Calculation showed [21] that results of this analysis are 
described by a deep attractive NN potential with one forbidden bound S-wave state. The 
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forbidden state has a wave function without a node. As a result, the 1 So-wave scattering 
wave function has a short-range nodal structure instead of short-range suppression specific 
to a repulsive co re p otential (RCP). After that, a preliminary attempt was made within the 
concept of MP [32] to reconstruct A^ A" potentials for the lowest partial waves (S and P) 
from data of the pp and pn PWA extended at the time to intermediate energies. 

At the same time the quark microscopic foundation of the MP remains the principal 
problem. Unlike the nuclear shell-model picture of the a — a interaction the lowest quark 
configuration s 6 is not forbidden by the Pauli principle and the corresponding OS-wave state 
of relative A^ A" motion is not forbidden either. Microscopic quark investigations of the last 
two decades with various kinds of qq interactions have resulted in the following short-range 
properties of the A^A" system 33J. There is a strong mixing of different six-quark configura- 
tions in the overlap region of two nucleons. For the S-wave states the leading configurations 
are s 6 and s A p 2 with comparable weights and destructive interference. This destructive inter- 
ference leads to strong short-range suppression of the A^ A" wave function. The suppression 
is described effectively by an RCP 34j. The s 4 p 2 [42] :r configuration introduced in our pa- 
pers [2lj and corresponding to the 2S-state of relative A^A" motion (i.e. to the MP) would 
dominate for instance in case of strong instanton induced quark-quark interaction but this 
interaction is not strong enough 35( . Further investigations 3g, |37| showed that there exists 



a source for strengthening of the s 4 p 2 configuration. Namely, if coupling of the NN, AA 
and hidden color CC channels is taken into account within the resonating group method 
then the symmetry structure of the highly dominant six-quark configuration s 4 p 2 implies 
the existence of a node in the S-wave relative motion wave function at short distances. Such 
nodes are specific to the MP. In the same manner microscopic qq interaction may give a 
short-range node in a P-wave of relative A^ A" motion wave function (in case of dominant 
six-quark configuration s 3 p 3 [33] x ). 

In summary, the question which type of the potential (MP or RCP) would be equivalent 
to the short-range quark microscopic picture of the A^A" interaction is highly controversial. 
For any RCP a phase equivalent supersymmetric partner with forbidden states (i.e. an MP) 
may be constructed |38|. Therefore these potentials are indistinguishable for the AW PWA. 
Specific to the MP appearance of short-range nodes in S- and P-wave relative motion wave 
functions is a result of complicated six-quark dynamics which is yet to be clarified. The nodal 
behavior of the MP wave function means that the wave function is not suppressed at short 



6 



range as in case of an RCP. Thus the MP produces high-momentum component richer than 
an RCP. This high-momentum component may be seen in electromagnetic reactions with two 
nucleons. In Ref. 39] it was shown that the available MP produces too rich high-momentum 
component in contradiction with the deuteron electromagnetic form factors. Thus we use 
the latest high-energy PWA data to refine short-range part of the MP. In our papers 20, |40 ] 
it was shown that the hard pp — > pp^ bremsstrahlung at moderate energies (Ei a f, ~ 500 
MeV) is critical to the kind of potential (MP versus RCP). The available experimental data 
at Sm a»er energy of ^ = 280 MeV fl g,ve only preliminary indication of MP validity 
[20J. Our present paper strengthens this line of phenomenological research using modern 
deuteron photodisintegration data. 



III. RELATIVISTIC OPTICAL NN POTENTIAL 



We apply the method of inversion |22| to the analysis of NN data up to energies at 
which relativistic effects are essential. We take into account these effects in the frames of 



the RQM Hi, ll5j. A system of two particles is described by the wave function, which is an 
eigenfunction of the mass operator M. In this case, we may represent this wave function 
as a product of the external and internal wave functions. The internal wave function \x) 
is also an eigenfunction of the mass operator and for system of two nucleons with masses 
mi = m2 = m satisfies the equation 

1 (1) 



M|x) 



2Jq 2 + m 2 + Vi 



int 



\X)=M\ X ), 



where Vi nt is an operator commuting with the full angular momentum operator and acting 
only through internal variables (spins and relative momentum), q is a momentum operator 
of one of the particles in the center of mass frame (relative momentum). Rearrangement of 



gives 



q 2 + mV 



X 



Q 2 X, 



where V acts like Vi nt only through internal variables and 

„ M 2 



Q 



- 2m 2 . 



(2) 



(3) 



Eq. d2J) is identical in form to the Schrodinger equation. The formally same equation may 
be deduced as a truncation of the quantum field dynamics [23J. The quasicoordinate repre- 
sentation corresponds to the realization q = V = V(r). 
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We applied the method of inversion 22j to reconstruction of the nucleon-nucleon partial 
potentials 

V(r) = (1 + ia)V {0) (r), (4) 

for single waves and 

v , r)= ( (I + ^OV/V) (l + ia 3 )V± \r)\ 

for coupled waves, where V^\r) are energy-independent real and inelasticity parameters a 
depend on energy. As input data for the reconstruction we used modern PWA data (single- 
energy solutions) up to 1200 MeV for isoscalar states and up to 3 GeV for isovector states 



of the NN system 



24j . The deuteron properties were taken from 25j. These data allow us 



to construct Moscow-type NN partial potentials sustaining forbidden bound states. These 
potentials describe part of the deuteron properties and the PWA data by the construction. 

According to the MP concept and the GLT some phase shift data of 24] are raised 180°. 
Namely, 1 So- wave phase shift and all four 25+1 P/-wave phase shifts are equal to 180° at zero 
energy, 3 Si-wave phase shift is equal to 360° at zero energy. The mixing parameters e% and 
e 2 of the MP differ from that of a traditional RCP by sign. All phase shifts for higher waves 
(for L > 2) are "small", they have zero values at zero energy. According to our model we 
have fitted free parameters of the inversion solutions to get nodes at r ~ 0.5 fm in S and P 
waves and to make central parts of the potentials close to each other and to the Gaussian 
shape. The energies of forbidden states are in the range 300 — 750 MeV. 

Our calculations show that the final state interaction (FSI) in the S and P waves gives 
by far the largest contribution to the deuteron photodisintegration cross-section comparing 
with FSI in other waves, so we present results of inversion only for these and for coupled to 
them waves. Part of results presented in Figs. [TH (for 1 Sq and 3 SD\ waves) we presented 



earlier in Ref. 



221. 



The reconstructed potentials V^°'(r) are displayed in Fig. [TJ The inelasticity multipliers 
a are displayed in Fig. [2j Fig. [3] displays reproduction of the corresponding phase shifts and 
mixing parameters. In Fig. H] the description of inelasticity parameters is shown. All the 
P-wave phase shifts are positive according to the GLT. Large difference between 3 Po-wave 
and 3 P2-wave phase shift curves reflects a large spin-orbital interaction which is attractive 
for the 3 P2-wave as we see. These features correspond to the general properties of the MP 
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FIG. 1: Reconstructed partial potentials for lower orbital momentum (single and coupled channels). 




FIG. 2: Reconstructed inelasticity multipliers a for the potentials presented in Fig. [TJ 

(its large positive gradient in the region r < 1 fm). It is interesting to learn from Fig. [3] that 
among four lowest pp phase shifts three of them ( 1 5'o, 3 Po an d 3 -P2) correspond to the MP 
but the experimental data within the energy range Ei a b = 2 — 3 GeV are contradictory for 
the 3 Pi-wave phase shift. It would be important to refine the PWA data in this range using 
modern polarization data on pp-scattering. The S- and .D-state wave functions of deuteron 
are displayed in Fig. [5j There is a node in the S'-wave function at r ~ 0.5 fm and both 
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wave functions are not suppressed at short-range in contrast with wave functions produced 
by an RCP. For continuum S- and P-wave functions the node radii equal to 0.5 — 0.9 fm at 
the considered energies. All potentials and inelasticity multiplies (a's) can be accessed via 
a link to the website 42]. 




?IG. 3: Phase shifts and mixing parameters in the present optical model. PWA data are from 



24]. For 1 Sq and 2S+1 Pj waves, the original data set from Ref. 24] is raised 180°. To leave the S 



matrix unchanged, we change the sign of the mixing parameters e\ and Si- 

It should be pointed out that in nuclear matter calculations the NN potentials should 
be used in the following form 

V nud (r) = V(°\r) + \(xs,L,j\ (6) 

where operator (xs,l,j\ projects onto the forbidden state \xs,l,j), positive constant A tends 
to infinity. The forbidden state \xs,l,j) m ay be found from Eq. (T5]) by some numerical 
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FIG. 3: (Continued). 

method as a bound state of the partial potential V^(r) (all bound states are forbidden 
except the deuteron one). Constant A is a large number, such that its further increase does 
not change results of calculation. This procedure orthogonalizes the nuclear wave function 
to forbidden two-nucleon states. Thus, we exclude the unphysical collapse of nuclear matter. 
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FIG. 4: Inelasticity parameters p in the present optical model. PWA data are from 24] 



IV. DEUTERON PHOTODISINTEGRATION IN POINT-FORM 
RELATIVISTIC QUANTUM MECHANICS 



Formalism of the PF is considered in detail in Refs. [Ill, [15[, while general covariant 
PF expressions for the electromagnetic current operator for composite systems are given in 
Refs. jlol . flf] ]. Therefore, we give only results necessary for our calculation, in notation of 
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FIG. 5: The deuteron S-wave and D-wave functions for present version of the NN Moscow poten- 
tial. 



Rcf. 



15] . We use algorithm of 15] for calculation of the matrix elements of the electromag- 



20| | . Similar approach 



netic current operator. We applied this formalism to th e py y process 
was applied to the elastic electron-deuteron scattering 19]. 

We consider the pn system and neglect difference of neutron and proton masses (mi = 
m 2 = m). Let pt be the 4-momentum of nucleon i, P = (P°, P) = pi + P2 be the system 
4- momentum, M be the system mass and G = P/M be the system 4- velocity. The wave 
function of two particles with 4-momentum P is expressed through a tensor product of 
external and internal parts 

\P, x ) = U 12 \P)®\x), (7) 
where the internal wave function \x) satisfies Eqs. ([TJ-tEl). The operator 

2 



U12 = Uu{G,q) = Y[D[Si;a(pi/m) 1 a(G)a(q i /m)} 



(8) 



i=l 



is the unitary operator from the "internal" Hilbert space to the Hilbert representation space 
of two-particle states [if] ]. _D[s;u] is the representation operator of the group SU(2) cor- 
responding to the element u GSU(2) for the representation with the generators s. Action 
of D[s;n] and matrices a are defined in Appendix A, Si — 1/2 is spin of a nucleon. The 
momenta of the particles in their cm. frame are 



qi = L[a(G)} x pi, 



(9) 
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where L[a(G)] is the Lorentz transformation to the frame moving with 4- velocity G 
(L[a(G)]~ l is the inverse transformation). It is easy to verify that qx = — q 2 . 
The external part of the wave function is defined as 

(G|P'>^^V(G-G'), (10) 

with scalar product 

(P"\P') = J ^(P"\G)(G\P') = 2y/M'* + P'*6 3 (P" - P'), (11) 

The internal part of the wave function \x) is characterized by momentum q = qi = — q2 of 
one of the particles in the cm. frame. Interaction appears according to the Bakamjian — 
Thomas procedure P = GM, where M is sum of the free mass operator M and of the 
interaction V: M = M + Vi n t (compare with Eq. ([I])). The interaction operator acts only 
through internal variables. Operators M, M, Vi nt and V commute with spin operator 5* 
(full angular momentum) and with 4-velocity operator G. Interaction term is present in all 
components of total 4-momentum. Generators of Lorentz boosts and generators of rotations 
are free of interaction. In the cm. frame the relative orbital angular momentum and spins 
are coupled together as in the non-relativistic case. Moreover most non-relativistic scattering 
theory formal results are valid for our case of two particles [11]. 



The deuteron wave function \Pi,Xi) is normalized as follows 

(p>,xi\p?,x*) = 2p? ^-pd- (12) 

For one-particle wave functions normalized in the same manner the free two-particle states 
are normalized as 

(p',x'\p",x") = ¥i\pWMW'i s m = 

= 4w(p' 1 )w(p' 2 )S 3 (p'( - pi)5 3 (p2 - p' 2 )<Vi/4' ( Ws ! ' = 
= 2W(P')6 3 (P" - PO^^W - ^8 M 8 M = 

= 2^(P')5 3 (P" - P')^5 3 (q" - d')8» W{ 8 M (13) 

where w(p) = ^m 1 + p 2 , M(q) = + q 2 , W(P) = VM 2 + P 2 , G (G) = y/l + G 2 , m 

are spin projections in the cm. frame. Multiplier is a relativistic invariant, therefore 
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we may normalize internal part of the scattering state wave function in the non-relativistic 
manner 



(-f ) X |-f i "X. )n.r. 

= 2^(P')o 3 (P" - P')5 3 (q" - q')<WW, (14) 

where S and \i are full angular momentum and its projection in the cm. frame. 
The differential cross section for the — > np process is given by 



dQ 64tt 2 M?/L 



Af\ 2 , (15) 



where g/ is the final asymptotic np relative momentum, k c is photon energy in cm. frame. 
The g?7 — > np amplitude Aif is defined in the same manner as the ppry amplitude that was 



used in Ref. 



20| 



(2ir) A 6 i (P t + k-P f )A lf = 
= V^jd A x (P f , X f\ e»Hx) \P h Xl ) e ikx (16) 

where Pi and Pf are initial and final 4-momenta of the NN system correspondingly, is 
the photon polarization vector. 

Following Ref. [ljj] we choose for calculation of the invariant amplitude a special 
frame defined by condition 

G, + G / = 0, (17) 

where G{ = Pi/Mi, Gf = Pf/Mf are 4-velocities of initial and final iViV cm. frames 
respectively (Gj and G/ are their 3- vector parts). The initial mass Mj is the deuteron 
mass. The final mass Mj is the invariant mass of the final NN system. These masses are 
different due to absorption of a photon, therefore coordinate frame corresponding to Eq. 
([17]) is not equivalent to the Breit frame where Pj + P f = 0. Masses Mj and Mf define also 
corresponding wave functions through Eq. ([3]) and Eq. (j2j). 

The matrix elements of the current operator J^(x) appears especially simple in the frame 
defined by Eq. (1T7I) : 

{P f , X f\Hx)\Pi,Xi) = 
= A^MMfe^f-^ ( X f\f(h) \ Xi )n,, (18) 
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where j M (h) is the current operator in the internal space defined in Ref. [15| as an operator 
J(0) (see Eq. (]A16j) ) in the frame ( 1TT|) expressed through h and q. Following Ref. [15J we 
use the dimensionless vector h = Gf/G®, where Gf is a 4- velocity of the final NN system 
in the frame defined by Eq. (fTTI) . This parameter may be expressed through the photon 
momentum k, so that h = 2(M i M f f/ 2 (M i + M f )- 2 k, |h| = h = (Mj- M f )/(M i + M f ) < 1. 
Convenience of this parameter is illustrated in Appendix B. 

The internal wave functions of the deuteron and of the final scattering state are normalized 
in the non-relativistic manner. The deuteron wave function is 

\Xi) = \Xi)n.r. = ~ E UjMMilMj), (19) 
r 1=0,2 

with normalization (xi\Xi)n.r. = 1; where 

\l, S; JMj) = E \S, /*> yUn) C^g, • (20) 
The internal wave function of the final continuum np state is 



2 1 °° 

\Xf) = \qf,3f,fif) n . r . = \/-— E x 

v ^ If j=o 

J J+5 J+S 1 I 

x E E E E »^/ff/.0O w 3^W/)|i , ,fi>;J r ^j>, (2i) 

Mj=-J l=J-S l'=J-S m=-l 

with normalization (x/'lx/)n.r. = <Kq} — Q/)^S/5' ^m/m' ■ ^ ne corresponding plane wave 
10/) n.r. is characterized by the spherical Bessel functions ji(qf,r)Sw instead of uf l ,(qf,r). 
The deuteron partial wave functions Ui(r) presented in Fig. [5] and partial waves of the final 
np states ufi,(qf,r) are calculated from Eq. (J2]). 
We define a reduced amplitude 

Tfi = (X/|<?(h) \Xi)n.r, (22) 



As a result the differential cross-section ( 1151) can be rewritten as 

da T^qfA^ 

where we average over photon polarizations, spin orientations of initial deuteron and sum 
over spin orientations of final nucleons. 
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In our calculations we approximate the above matrix element 



( X f\ f(h) \ X i)n,. « (0/1 J^h) \ X i)n,. + (Xf - 0/1 J (h) |Xt)n.r. • (24) 

The first term is a plane wave approximation (P1WA) and is calculated using the exact 
current operator (1A6I) . In this case the operator q can be substituted by q/ and operator 
structure of j^(h) can be presented as 

Hh) = /(h) + 5f = £ (Bfi + (B£ • Si ) + (B& • s fc )+ 



i=l,2 



+ (B^-s i )(B^-s fc ))/ l (h)+5j, (25) 

where j M (h) is sum of the one-nucleon electromagnetic current operators (spectator ap- 
proximation), addend (5j M restores the current conservation equation; k = 2, if i =1, and, 
conversely, k — 1, if % — 2. and B^j , m > 2 are vector and tensor functions of arguments 
h and q^. These functions are given in Appendix A. In Appendix A we calculate addend 
if from the current enervation sanation fC.owing Ref. Q as we did for the W , proeess 
in Ref. [20]. Obviously our phenomenological quasipotential model offers no microscopic 
picture of interaction that would allow us to unambiguously determine the current opera- 
tor. We use the defined bellow <5j M only to estimate violation of the current conservation 
equation. Assuming gauge invariance (which follows from the Poincare invariance and the 
current conservation equation) we use the transverse gauge 

e M =(0,e), (ek) = 0. (26) 

Thus, we exclude the j°(h) and j||(h) (see Eq. (1A21D ) components of the current from 

)y definition of the current operator J M (a;) 



m) 



Eq. (|22|1 . The Poincare invariance is ensured 
through the operator j^{h) (see details in Ref. 

The use of (x( r ) — 0( r )) combination in Eq. ([21]) accelerates convergence of the partial 
wave expansion. This term is nonzero due to FSI of neutron and proton. It is calculated 
from the first order in h approximation of the current operator -Wh). This approximation 



calculated in the same manner as one for the pp system in Ref. 20J is given by 



j(h) « j(h) = 6j + ^f(O) - hGT(O) + 

+i f-[S x h] + — -[q x h] (q ■ S)) GZ l (0) + 

\w L 1 w{w + m) 1 J v ' j mV 1 
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m 



+1 -[Txh] 



+ m) 



[q X h] (q ■ T) gf-'(O) + 



l(h . [qxS])q ^ + 6r(oi 



.(h.hxTDqf^ + ^L' 

y mw w{w + m) 



\ raw w(w + m) 
GT(0) 



+ 



(h • q) q 



(27) 



<J = f TT~T~7T" - 1 - h) V«» + 

\ Mf + Mi J w 



+lh [q x T] 



/GT(0) GT(0)\ 



V m w + m J 
where S = Si + s 2 , T = s x — s 2 , 



2f(0)«;r+[qxS] 



g?{0) = Gl (0)/i(h) - G?(0) J 2 (h) , 
C n (0)=Gf n (0)/ 1 (h)-G^(0)/ 2 (h), 
G'n0)=^(0)/ 1 (h)+G^(0)/ 2 (h), 

Gr(0) = G^(0)/ 1 (h)+G:(0)/ 2 (h) 



u> = to(q) 



m 2 + q 2 , 



u> + m 



(28) 



(29) 
(30) 



di(q) = q+(-i; 



2h 



1 - /i 2 



/i(h) X (q) = X(di(q)), 
(tw + (-l)^h ■q))«q+ (-l) < 2hw . 



(31) 



where G^Qjy), G^{Q 2 N ), G™(Q%), G p e {Q 2 N ) are nucleon electromagnetic form factors pa- 



43) 



rameterized according to Ref. 

In Ref. 1q| the elastic electron-deuteron scattering was described in frames of the PF 
RKM. It was shown that in the PF SA that the momentum of the unstruck particle (the 
spectator) is unchanged, while the impulse given to the struck particle is not the impulse 
given to the deuteron. 

Following a general approach to construction of the electromagnetic current operator for 



relativistic composite system 15] we define t 
an increment of the particle 4-momentum q$ 



re momentum transfer Q\ to the particle % as 



19] 



Q 2 i = M - Qi 



(32) 
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For interacting particles the individual 4-momenta are not defined before photon absorption 
as well as after it. Therefore we introduce an operator Q 2 corresponding to the physical 
quantity of the momentum transfer Q 2 . In Appendix B we generalize the deduction presented 
in Ref. [3] and show that 

Qt = - ?i) 2 = 16(- 2 + q 2 - ^!) JT ^y r (33) 

This is the general expression of the Q\ = Q\ = Q 2 N in case of free two-particle states (for 
particles of equal masses), therefore we use this expression in the PF SA for evaluation of 
the current operator in Eq. (I25j) . The parameter h does not depend on interaction and is 
specified by relative "position" of initial NN cm. frame and final NN cm. frame. In 
case of two interacting particles q and Qf are operators in the internal space. In impulse 



representation q is a variable of integration 19|. It is obvious that in action on a plane wave 
(for P1WA) this operator is equivalent to the multiplication by a number Q 2 N > if h ^ 0. 
The first order in h approximation gives Q 2 N ~ for Eq. ( 1271) . Consideration similar to that 
of Ref. Q gives for the P1WA in our case of the deuteron photodisintegration 

Q 2 N = E 2 - (w' n - w' p f = (2w' n - m D )(2w' p - m D ), (34) 

where w' n and w' n are final energies of neutron and proton in the initial cm. frame (laboratory 
frame). Detailed deduction of Eq. (1341) is published in Ref. [44]. 



V. RESULTS AND FUTURE PROSPECTS 

Our theoretical description of the differential cross-section of dj — > pn reaction is com- 
pared with recent experiment jj] in Figs. [MS at a few energies around E 1 = 2 GeV. We do 
not use free parameters. However, there are uncertainties in our calculation. First uncer- 



tainty is caused by uncertainty in the form factor parametrization of Ref. 43j due to errors 



of the experimental data on the form factors. We estimated this uncertainty at about ±15 



per cent of the results varying parameters inside limits defined in Ref. 43j. Second uncer- 
tainty is connected with the approximation (1271) used to calculate the FSI term in Eq. (|24l) . 
In Eq. (|27p the nucleon form factors are equal to their values at Q 2 N = and therefore the 
FSI term is overestimated. Fig. [7] (J5 7 = 2.051 GeV) shows contribution of the FSI term. 
The P1WA term of the amplitude (|2"4"|) is dominant, but the FSI is not negligible. Therefore 
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it is desirable to estimate the second order in h correction to the approximation ([27]) . We 
plan to do this estimation in the future. Third uncertainty is caused by uncertainty of the 
addend <5j that restores the current conservation equation. To estimate this uncertainty 
we calculated two curves for every energy of the photon. The lower curves correspond to 
calculations without addend <5j in the current operator and with form factors of Ref. 43 ] 
varied to their lower limits. The upper curves correspond to full calculations and with form 
factors varied to their upper limits. The FSI is included for both curves. 

We see good general correspondence of the theory and experiment both in absolute values 
and in shape of angular dependence of the differential cross-section at various energies. Large 
absolute values of cross-sections in our theory in comparison with results for the RCP (Fig. [7J 
(Ej = 2.051 GeV)) originate mainly in the nodal character of the deuteron 5*-wave functions 
(greater weight of the high-momentum wave function components). The ability to describe 
both the absolute value and the angular dependence of differential cross-sections confirms 
the detailed algebraic structure of our theory. A persistent forward-backward asymmetry is 
determined mainly by the angular dependence of the nucleon electromagnetic form factors 
according to Fig. [9] (proton knockout dominates at forward angles and neutron knockout 
dominates at forward angles). 

To complete this line of our investigation we plan to make an analysis of polarization 
d-1 ~> P n experiments and to consider the pionic radiative capture pp — ► d-K + at proper 



energies. Other actual problems are outlined in Ref. 20]. First of them concerns the 
microscopic theory of the MP. As we suppose, it is connected to the short-range quark 
exchange between nucleons accompanied by excitations of color dipole states of two virtual 
baryons with very strong attraction between them. This scenario is based on the quark 
configuration s 4 p 2 [42]a. in deuteron. 

As the concluding remark it should be stressed that usage of the MP instead of an 
RCP in the theory of complex nuclei demands accurate evaluation of 3N— forces. Effect of 
these forces is much enhanced j^j], as far as three nucleons without NN core can overlap 
and form short-range 9g-subsystems with large probability. Recent experiments [46( on the 
knock-out of nucleon from 3 He nucleus may clarify the situation. In these experiments the 
missing momentum is great and recoil to 2A-subsystem with large relative momentum of 
two spectator-nucleons is observed. 

Authors are grateful to Prof. V.I. Kukulin and to Dr. N.A. Burkova for useful remarks. 
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O 30 60 90 120 150 180 O 30 60 90 120 150 180 



e p scm (degrees) 6 p scm (degrees) 



FIG. 6: Angular dependence of d'j — > pn reaction differential cross-sections for different photon 
energies E~. Our theory is compared with experiment 
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FIG. 7: The same as in Fig. [6j but (for E 1 = 2.051 GeV) our theoretical results for the RCP (Paris 
potential [45J) is also shown (two lower dashed curves). Upper dashed curves show results of our 
calculations for P1WA with deuteron wave function in the initial state calculated with our MP. 
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FIG. 8: The same as in Fig. [6] 




FIG. 9: Angular dependence of 4-momentum transfer Q 2 , and of the nucleon electromagnetic form 
factors for d'y — * pn reaction calculated from Eq. (I34h for the P1WA. In our calculations we use 
dependance of the form factors on Q 2 according to parametria of J.J. Kelly Q. 

APPENDIX A: 



In this Sec. we explain calculation of the electromagnetic current matrix elements. The 
derivation is based on results of Ref . [15j , where Eq. (ITS]) and Eq. (1A6I) were deduced. 
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Let us dolino a. matrix Ji 

= / = = Al 
12 + l 



corresponding to a 4- velocity g>, where cr = (cr^, a y ,a z ) are the Pauli matrices. Let us define 
the matrix p = M(p) = a^p^ corresponding to a 4- vector p (a is 2 x 2 unit matrix). 
Operator M(p) transforms the 4-vector p to (2 x 2) matrix. The inverse transformation is 
defined as 

P0= ^{Pn+P22), Pi = „(Pl2+P2l), 

P2 = ^t(P21 -Pl2), P3 = g(Pll -P22), (A2) 

we denote this transformation as p = V(p). The boost p — > is equivalent to matrix 

transformation 

p a(g)P a (g) + - (A3) 

It is easy to see that L[a(g)](l, 0, 0, 0) = g. The Poincare group transformation U(a,l) is 
characterized 15( by the 4-shift a and 4-rotation 

U(a, l)ip{g) = e m °' a D[s; aig^lafflltptf), (A4) 

where <f{g) is a normalized spinor function of a particle with mass p; s is spin of the particle; 
and g' = L(l)~ 1 g. In our case of spin s = 1/2 particles, we deal with the fundamental 
representation [43], i.e. s.; = ~<Ji and 

D(s; a(g)-Ha(g')) = a(g)- l la(g') . (A5) 

The "internal" electromagnetic current operator for a system of two particles in the SA 
is m 



where 



j»(h) = ^ (V)*D[Dtir{h)DiK i Ii(h), (A6) 

i=l,2 



(^^l^W/)]!:,^/')]^:)", (A7) 



D\ = D\s k \a(q k /m k ) 1 a{f) 1 a(f')a(d ki /m k ) 

= a k {q k /m k y 1 a k {f)'' 1 a k {f)a h (d ki /'m k ), (A8) 
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D\ = D 



Si]a(qi/mi) l a(f) 1 a(z i 
= a i {q i /m i )~ 1 a i {f)~ l a i {z i ) , (A9) 



Dl = D 



Sj 1 a(z i ) 1 a{f)a{di/mi) = 

= ouifi)- 1 ^ {ziT 1 on{f)ai{di/mi), (A10) 



kinematical multipliers 

3/2 



r =wSJ (A11) 

Here, k = 2, if z =1, and, conversely, = 1, if % — 2. L{G, G') denotes the Lorentz transfor- 
mation L [a (G, G% anda(G, G') = a((G + G')/\G + G'\), z { = L[a(/)]g i /m i , L[a(f')]di/mi. 
Next, 

/ = L(G, G') _1 G, f = L(G, (A12) 

represent the 4-velocities of the two-nucleon cm. in the initial and final state, respectively, 
meaning the coordinate frame f|T7|) . The following formal aspects should be mentioned here: 
P = f 2 = 1, f+f = 0,f° = f = (1+f 2 ) 1 / 2 , h = f//°; L(G, G') = L(a(G, G% a(G, G') = 
a((G + G')/\G + G'\); d x = ( Wl (di),di), d 2 = (w 2 (d 2 ),d 2 ), d 12 = L[a(f )" 1 a(/)]g 2 = 
(ty 2 (di), — di), d 2 \ = L[a(f')~ l a(f)}qi = (wi(d 2 ),d 2 ); last equations gives dj also. Index z 
or k of matrices a and a means that it acts in z-th or k-th particle spin space and appears 
as in Eq. (1A1I) but with cr, or au correspondingly instead of a. Let d\ = (a>i(di), di), c? 2 = 
(u; 2 (d 2 ), — d 2 ) and Ij(h) (z = 1, 2) be operators defined by the conditions ij(h)x(q) = x(cli)- 



;q.) /M(d 



Finally, .77(h) is a 4-current of the particle i, 



j°(h) = eF* (Q 



2 X 



m = ~^—F l m (Ql) x Si ), (A13) 



where vectors hj are defined below, Sj = Cj/2, Q 2 = 4m 2 h 2 /yl — h 2 (see also Eq. (1551) ). 

From (1A6HA131) it is obvious that for a plane wave final state, when operator q = — z'V can 
be substituted by vector q^, operator j(h) becomes an exterior product j u (h) = J2i=i 2 
Aj^I^h); fc = 2, if z =1, and, conversely, = 1, if z = 2. The q/-dependent matrix A^j, 
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acts in i-th particle spin space and presentation A k v = (y^or ivVi = a k uQ + 2(sj • & k u ) is valid. 
"Components" a u are extracted by (IA2[) transformation 

aU = V((L%D^(h)DlK>) 

a k ui =V(Dl), i^k. (A14) 

Functions B of Eq. ( 1251) are expressed as 

Df n i n k T3V n k _i of n i k 

k = 2, if i =1, and, conversely, k = 1, if i = 2. 

Now, we should take into account the current conservation equation 

Using also the 4-shift 

= exp{iPx)J>*{0)exp{-iPx), (A16) 

we obtain a relation 

P M >(0) - >(0)P M = 0. (A17) 

In terms of the internal variables of iViV-system, Eq. ( 1A17I) can be reduced to the matrix 
element 

< X f\M f G fo f(h) - M f G f j(h) - MiG io f(h) + MiGa(h)\xi >= 0, (A18) 

that can be rewritten in the form 

< X /|(h • j(h))| X< >= < XflhVlXi >, (A19) 

as far as G; = -hG l0 , G, = hG i0 , P { = MjGj, P, = M f G f , M\ Xi >= Mi\ Xi >, M\Xf >= 
Mf Wf >■ The current (1A6j) does not satisfy Eq. (1A22[) and needs a modification. Following 
15j we use the unique decomposition into longitudinal and transverse parts: 

h 

h' 



j(h)=j(0) + -j||(h)+j ± (h), (A20) 



where hjj_(h) = and 

J||(h) = ^ T (h-(j(h)-j(0))), 



j ± (h) = j(h) - j(0) - A (h- (j(h) - j(0))) . (A21) 
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To estimate violation of the current conservation equation we assume that NN interaction 
does not change transverse and time components of operator j(h). Then we can reconstruct 
j(0) and j||(h) from Eq. ( 1A22I) . In the transverse gauge (1261) the longitudinal component 
has no effect on results of our calculation and therefore we determine only matrix element 
ofj(O) 

\Xi >, (A22) 



Mf-Mi . 8j°(h) 

Corresponding addend <5j that restores Eq. (1A22p is given in Eq. 
The first term in Eq. ([21]) (P1WA) appears as 



(0/l?lxi>».r. = v~i; E E i l (^(i2j2yL(W,s;JMj\ii\i,i- } iMj)ut 

V 71 Qf J=o 1=0,2 m=-l \k=li=l 

+3U?/) (1 - <W) (£<J, S; JMj\K?\l, 1; lM,)^) 
\j=i 

/■OO 

-2g?w(q f ) "£(l',l;JMj\f%l',lMj) r>( g/ r>,(r)efr 

l'=l,3 J ° 

^M?/) = / ji(qfr)ui(r)dr, 



o 



^ - mfa/)(1 ; ( g) )( +(-i)'2(h. q/ ) r^w^M*- < a23 » 

here Ljj^ (A; = 1, 2, 3, 4) = B^, (B^s 2 ), (B^Si) and (B^^XB^Si), respectively; K? represent 
the /i-components (p = 1, 2, 3) of the the first three (i = 1, 2, 3) terms in Eq. ([28]) . 
The second term in Eq. (I24p appears as 



oo J+S J+S I 

(xf-<t>f\T\K)\xi)n, = \l~ E E E E E^PU?/) 

"5/ L=0,2 J=Ql=J-Sl'=J-Sm=-l 



r°° , a 

x ^ dr(Z',S; JM|(^( g/) r)-Vjj(?/r))j (hK(r)|L, 1; 1M,) (A24) 

We use further the algebraic results (A24)-(A27) of Ref. [^(J and obtain the final ex- 
pression for the differential cross-section which is reduced to radial integrals and spherical 
harmonics but, unfortunately, is too much unwieldy to be exposed here. 

Now, by a few examples, we illustrate the calculation technique for the matrix elements 
of various components of the relativistic current operator 

(l f ,S f ;J f M f \(B^- Sl )\l h S i ;J i M i ) = 
= (B$ i .(l f ,Sf,JM f \a 1 \l i ,S i ',J i M i )), 



(If, S f ; JM f \(si) v \li, S^, J t Mi) = 

= (- 1 ) Li+Jf+Sf+ls L i L f C(i / ^ fMf ^2J i + 1 x 
x \L f J f S f L 5 ||aill5i)j 

(SfWsiWSi) = (-1)^(2^ + l)(6S / + 3)/2x 

x < 1/2 1/2 S f 
1 $ 1/2 



(ai ■ Si)(a 2 • s 2 ) 



k=0 



J2 C k [a x x a 2 ] ik) x [si x s 2 



(o) 



C k = (1, v^, >/5); 



(5 / ||[ ai x a2 ](*)||5 i > = 



1^/(25* + l)(25 / + l)(2fc + l) 



1/2 1/2 5/ 
1/2 1/2 5, 
1 1 k 



(V • S)V A 



V3 



[V x V] (0) x 5 



(i) 



[V x V] (2) x S 



(i) 



—2- 



(h-[VxS])V M 
V6 /-s/15 



[[V x V] (2) x Sp) x /i 



(i) 



+ 



[[V x V] (0) x x h 
^[[[VxVfx^Wx^ 1 ' 



(i) 



•, S f = 1; J / M / | [V x V]<*> x S]« W /(r) | L,, S t = 1; J,M 4 > 



1 ■// 



Lj 1 Jj 



kin 



y / 6(2J l + l)(2n + l) < L^lfV x Vp/Ml^i >, 
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<L f \\[Vx Vp^||L,> 



2L f + 1 l L , 3(2L 2 + 2L,- 1)^/2(2^ + 1) \ /rf 2 Li{Li + i)\ 

- S LiLf -1 + — lWor ; * r ; , ho 3 /( r ) 



v^Cf'° r \ ^ \ (2L t - 1)(2L, + 1)(2L, + 3) I 



3(^ + 1)^ + 2)^2(2^ + 1) /^ (2Lt + 3) rf (L i + 3)(L< + 1) \ 
+0L i% -2 (2L i + l)(2L i + 3)(2L i + 5) Ur 2 r dr r 2 J /lrj 



, r 3^-1)^2(2^ + 1) ( d * (2L.-1) d L^-2) ^ 

+< W 2 (2L, + l)(2L i - 3)(2L, - 1) ^ ~ rf7 + J /(r) ' " (A30) 

In these expressions, an upper index in round brackets means a tensor rank of an operator. 
The first rank is omitted where it is obvious (V = etc.) 



APPENDIX B: POINT-FORM MOMENTUM TRANSFER 

In the general case there are initial iViV-state with associated initial cm. frame (i.c.m.f.) 
and final iVTV-state with associated final cm. frame (f.c.m.f.) Suppose that the photon 
momentum (momentum transfer) is along the z axis. Values of photon momentum and 
energy in i.c.m.f. are |q 7 | and q® correspondingly. Momentum transfer is Q 2 = |q 7 | 2 — (f^) 2 - 
Let P be the total 4-momentum of the iViV-system, M be the mass of the iVTV-system, 
G = P/M be the system 4-velocity. Index i(f) means initial (final) state of the iViV-system. 



Transformation from i.c.m.f. to the special frame suggested by Lev 15] (L.s.) where 

G f + Gi = 0\ L . s . (Bl) 

is defined by angle A/2 such that 

tanhA/2 = /i, (B2) 

where h = G//G*|l. s .. The Lev frame ( IBlj) is not equivalent to the Breit frame defined by 
the condition Pf + Pj = if Mf ^ Mf. In case of elastic electron-deuteron scattering these 
frames coincide. 

From this point we may use a special derivation of Ref. [3] (Eqs. (IB3HB7I) of the present 
paper). 

The initial energies and z-components of momenta in L.s. are 



W\ = w cosh A/2 + q z sinh A/2 
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Qiz — Qz cosh A/2 + it; sinh A/2 
w 2 = w cosh A/2 — q z sinh A/2 
q 2z = - q cosh A/2 + w sinh A/2, (B3) 



where q and w = y/q 2 + m 2 are center of momentum variables, q is momentum of particle 
one (internal variable). After the photon absorption the z-component of the internal variable 
and corresponding energy change 

q' z = q z cosh A =F w sinh A (B4) 
w' = w cosh A =f q z sinh A, (B5) 

where the minus (plus) sign is used when particle one (two) is struck. The final energies and 
momenta in L.s. will then be 

w[ = w cosh 3 A/2 — q z sinh 3 A/2 
q' lz = q z 3 cosh A/2 — wsinh3A/2 

w' 2 = w 2 ; q' 2z = q 2z , (B6) 

other components do not change. Some hyperbolic trigonometry reveals that 

(q[ - qi ) 2 = A(q 2 z - w 2 ) sinh 2 A, (B7) 

it follows from Eq. (|B2]) that 

2h 

sinhA = - -. (B8) 

1 — h A 



Since 



the resulting Eq. f[33|) is established 



ql ~w 2 = -(m 2 + qi) = -{m 2 + q 2 - ^-), (B9) 
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